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Abstract. A very general surface of degree at least four in P 3 contains no curves other than 
Q ^ intersections with surfaces. We find a formula for the degree of the locus of surfaces in P 3 of 

■ degree at least five which contain some elliptic quartic curve. We also compute the degree of the 

Xf*^ ' locus of quartic surfaces containing an elliptic quartic curve, a case not covered by that formula. 



1. Introduction 

The Noether-Lefschetz theorem asserts that all curves contained in a very general surface F 
of degree at least four in P 3 are complete intersections. This is usually rephrased saying that 
the Picard group is Z. Noether-Lefschetz theory shows that, roughly speaking, each additional 
generator for PicF decreases the dimension of the locus of such F in F N = \O p3 (d)\, d > 4. 

Let W be a closed, irreducible subvariety of the Hilbert scheme of curves in P 3 with Hilbert 
polynomial Pw(t)- Let us denote by NL(W,d) the subset of P^ defined by the requirement that 
the surface contain some member of W. 

The purpose of this note is to address the question of determining the degree of NL(W, d) for 
the family of elliptic quartic curves in P 3 . 

When W is the family of lines, or conies, or twisted cubics, formulas for NL(W, d) have been 
found in [9]. There as here, we follow the strategy of using Bott's formula as explained in [5]. We 
get a polynomial formula (4.3) valid for d > 5. We also compute the degree (38475) of the locus 
of quartic surfaces containing an elliptic quartic curve. The case of quartic surfaces is not covered 
by the formula essentially because the map that forgets the curve shrinks dimensions: gcncrically, 
it contracts a pair of disjoint pencils, see 4.2. 

2. THERE IS A POLYNOMIAL FORMULA 

Let W be a closed, irreducible subvariety of the Hilbert scheme of curves in P 3 with Hilbert 
polynomial Pw(t). Let 

(2.1) WxP N P2 y ¥ N 

Pi 
W 
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denote the projection maps from W x P N . Castelnuovo-Mumford regularity [12] shows that for all 
d >> 0, the subset NL(W, d) of pairs (C, F) in W x P N such that the curve C is contained in the 
surface F is a projective bundle over W via p\. We have 

codim WxP nNL(W, d) =pw(d). 

For instance, if W is the Grassmannian of lines in P 3 , then pw(d) = d+1 and so dim N L(W, d) = 
N- (d-3), d> 1. 

Let us denote by NL(W, d) the subset of P N defined by the requirement that the surface contain 
some member of W. In other words, with notation as in (2.1), 

NL(W, d) =p 2 (NL(W,d)). 

We assume henceforth that the general member of W is a smooth curve. 

2.1. Proposition. For fixed W we have that degNL(W,d) is a polynomial in d of degre < 
3 dim W, for all d » 0. 

Proof. Let C C W x P 3 be the universal curve. Likewise, let F C P^ x P 3 be the universal 
surface of degree d. Write C,F for their pullbacks to W X ¥ N x P 3 . We have the diagram of 
sheaves over Y := W x P N x P 3 , 

O y ^O y (F) ^O p (F) 




5 (F) 

By construction, the slant arrow p vanishes at a point (C, F,x) £ W x P^ x P 3 if and only if 
x e F n C. We have G C F when the previous condition holds for all x € C (point with values 
in any C-algebra). Thus NL(W, d) is equal to the scheme of zeros of p along the fibers of the 
projection p 12 : C -> W x P N . Recalling [1 , (2.1),p. 14], this is the same as the zeros of the adjoint 
section of the direct image vector bundle p\2*{0 'q(F)) . Let 

W x P 3 ^ P 3 

gi 

w 

denote the projection maps from W x P 3 . Since O(F) = O p jv(1) <S) Op3(d), by projection formula 
we have to make and do with a section of Opjv(l) £g> £d, where 

(2.2) S d = qi *{O d {d)). 

By Castelnuovo-Mumford and base change theory, there is an integer do such that £d is a vector 
bundle of rank pw(d) for all d > do (^regularity, see Remark 2.2.). In fact, it fits into the exact 
sequence of vector bundles over W, 

(2.3) ^ qi *{l d ®O v3 {d)) »-gi*(Ppa(d)) ^<H*(O d (d)) ^0 

T^d > Td »- £d 

where we set for short 

F d = H°(O p3 (d)), 



SURFACES CONTAINING AN ELLIPTIC QUARTIC CURVE 



3 



(trivial vector bundle with fiber) the space of polynomials of degree d. Taking the projectivization, 
and pulling back to P N x W, we get 

<V(-i) 



T^d > Fd *»- £d 

By construction, p vanishes precisely over NL(W, d). This shows that we actually get 
(2.4) NL(W,d) =F(V d ). 

Since rank of £d and codimension of NL(W, d) agree, it follows that NL(W, d) represents the top 
Chcrn class of O p n (1) ® Ed (cf. [6, 3.2.16, p. 61]). This is the key to the calculation of degrees below. 
The map 

fot .s W x¥ N D NL{W, d) NL(W, d) C ¥ N 

[2 - b) (C,F) i — ^ F 

is generically injective by Nocther-Lefschetz theory [ ], cf. Corollary 4.2 below. Therefore the degree 
of NL(W, d) can be computed upstairs. Namely, setting 

m = dimNL(W,d),H = Cl O rN {l), 

we have 



degNL(W,d) = j H m n NL(W, d) = J H m c T (0 P N(l) ® £ d ), 
where r ^rankf^. Expanding the top Chcrn class and pushing forward to W, we arrive at 



(2.6) deg NL(W,d)= / c w (£ d ), 

Jw 

with w = dim W. Since £d is the pushforward of a sheaf on W x P 3 , we may apply Grothendieck- 
Riemann-Roch [6, p. 286] to express the Chern character of £d as 

ch{£ d ) = ch(q v .(O d (d)) = q u (ch(0g) ch(C» P 3(d))toddP 3 ) . 

Notice that the right hand side is a polynomial in d of degree < 3. Since c w is a polynomial of 
degree w on the coefficients of the Chern character, we deduce that c w (£d) is a polynomial in d of 
degree < 3w. □ 



2.2. Remarks. (1) The assertion that £d, as defined in (2.2), is a vector bundle of rank pw{d) 
holds for all d beyond do = the maximal Castelnuovo-Mumford regularity of the members of W. 
For instance, if W is the family of lines in P 3 , then do = 1. 

(2) For the case of elliptic quartic curves presented below, we note that the regularity of the ideals 
(x\ : x\) and {x\X2, x\, x\) is 3, whereas for (xq, XqXi, Xqx\, x\) it is 4. The last two ideals are 
representatives of the closed orbits in W . An argument of semi-continuity shows that d = 4 works 
for all members of W , see [2]. Nevertheless, the map in display (2.5) is generically injective only 
for d > 5, cf. (4.2) below. Notice that the full Hilb 4t P 3 has a "ghost" component with regularity 
index > 5, see [7]. 
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3. ELLIPTIC QUARTICS 



We consider now the case of surfaces of degree > 4 containing an elliptic quartic curve in P 3 . 
Thus, a general member C4 of W is the intersection of two quadric surfaces. The parameter space 
W is described in [2] and has been used in [5] to enumerate curves in Calabi-Yau 3-folds. For the 
convenience of the reader, we summarize below its main features. 

The Noether-Lefschetz locus of quartic surfaces containing some C4 is slightly exceptional. This 
is a case when the map (2.5) fails to be generically injective (cf. Corollary 4.2): it actually shrinks 
dim NL(W, d) = 34 to dim NL(W, d) = 33. Indeed, if a quartic surface F contains some general 
elliptic quartic C4, then F must contain the two pencils \C±\ and \C' 4 \, where C' 4 is the residual 
intersection of F with a quadric containing C4, i.e., 2H = C4 + C' 4 , with H = plane section. We 
show in §3.2 that NL(W,4) is a hypersurface of degree 38475 in P 34 = |Ops(4)[. 

3.1. Next we give an outline of the calculation. Put 

X = G(2,J- 2 ), 

the Grassmannian of pencils of quadrics in P 3 . 

The diagram below summarizes the construction of W as explained in [2] . 



(3.1) 



G(19, J4) 



^3) x 



D 



w 



D 



<G(2,.F 2 ) 



X 



X 



D 



D 



E 



D 



Z D 



Y 



Y 



where 



Z = P 3 x G(2, F\) consists of pencils with a fixed plane; 

Y = {(p,£) \ pD £} = closed orbit of Z; 

Y — > Y = P 2 — bundle of degree 2 divisors on the varying 
X = the blow-up of X along Y and 

X = the blow-up of X along Z. 



z 



D Y 



Y 




Let 

AcF 2 xX 

be the tautological subbundle of rank 2 over our Grassmannian of pencils of quadrics. There is a 
natural map of vector bundles over X induced by multiplication, 

H3 ■ A <g> F\ — > F 3 x X 

with generic rank 8. It drops rank precisely over Z. It induces a rational map k : X ---> G(8, J-~s). 
Blowing up X along Z, we find the closure X C G(8, F3) x X of the graph of K. Similarly, up on X 
we have a subbundle B C F3 x X of rank 8 and a multiplication map 



Hi : B <S> Fi 



F 4 xX 
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19 ~ 

with generic rank 19. The scheme of zeros of /\ \i\ is equal to Y . Indeed, it can be checked that 
each fiber of B is a linear system of cubics which 

• either has base locus equal to a curve with Hilbert polynomial p w (t) = 4i 

• or is of the form p ■ *, meaning the linear system with fixed componentc a plane p, 
and J~2* denoting an 8-dimensional space of quadrics which define a subscheme of p of 
dimension and degree 2. 

The exceptional divisor E is a P 8 -bundle over Y. The fiber of E over (p, 2/1+1/2) € Y is the system 
of quartic curves in the plane p which are singular at the doublet yi + 1/2 ■ Precisely, if xo, . . . , X3 
denote homogeneous coordinates on P 3 , assuming p := xa,t — (xq.xi), a typical doublet has 
homogeneous ideal of the form (xq, Xi, f(x2, X3)}, with degf = 2. Our system of plane quartics 
lies in the ideal (x\,f) = {x\, Xif, f 2 ). Given a non-zero quartic g in this ideal, we may form 
the ideal J = (xQ,XoXi,xof,g). It can be checked that the J contains precisely 19 independent 
quartics and its Hilbert polynomial is correct. In fact, any such ideal is 4-rcgular (in the sense of 
Castelnuovo-Mumford). Moreover, up on X we get a subbundlc 

of rank 19. Each of its fibers over X is a system of quartics which cut out a curve with the correct 
Hilbert polynomial. The multiplication map 

C <X> F d -A — ► Fd x X 
is of constant rank ( d ~^ 3 ) — 4d. The image 

is a subbundle as in (2.3). We have 

NL(W, d) = F(V d ) cf xi 

Now the map NL(W 1 d) — > NL(W,d) is generically injective for d > 5 in view of Corollary 4.2 (ii) 
below. Hence the degree of the image NL(W, d) C P N is given by J deSd, cf. (2.6). 

The above description suffices to feed in Bott's localization formula with all required data. 
Indeed, X inherits a C*-action, with (a lot) of isolated fix points. The vector bundle T>d — > X is 
equi variant; ditto for Ed- Bott's formula reads [4], 



(3.2) c w £ d = C -^f^- 

pGnxpts ±u ^ 

The equivariant classes on the r.h.s are calculated in two steps. We set below d = 5 for 
simplicity. First, find the C*-weight decomposition of the fibers of the vector bundles Ed and T at 
each fixed point. Then, since the fix points are isolated, the equivariant Chern classes cf are just 
the symmetric functions of the weights. 

For instance, for the tangent bundle T, say at the fixed point corresponding to the pencil p = 
(xq,x 2 ) g X, we have 

T p = Hom(p, Tijv) = p y <g> Tijv = ^ + ^ + ^ + -- - + 4- 

x Q x Q x x 1 

Each of the 16 fractions , f^j, on the right hand side symbolizes a l-dimensional subspace of 
T p where C* acts with character t x i+ x i~ x k-xi _ rp ne denominator in (3.2) is the corresponding 
equivariant top Chern class, to wit 

(xi - x )(x 2 - X )(X! +x 2 - 2x ) ■ ■ ■ (2x3 - 2xi). 
The numerator in (3.2) requires finding the monomials of degree 5 that survive modulo the ideal 



(i 
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(x%,xf) ■ (x , ■ • . ,a;3) 3 - 

We are left with 20 (=rank^5) terms, 

x 2 XiXq + X3X2X1 + • ■ ■ + xfxo- 
Now the cquivariant Chern class cj 6 (£d)p is the coefficient of t 16 in the product 
(t + xo + xi+ 3x 2 )(t + x\ + 3x 2 + x 3 ) ■ ■ ■ (t + x + 4x 3 ) 
(20 factors.) In practice, all these calculations are made substituting Xi for suitable numerical val- 
ues, cf. the computer algebra scripts in [13]. 

3.2. the case d = 4. Presently p 2 ■ iVX(W, 4) — > NL(W, 4) is no longer generically injective. It 
shrinks dimension by one: a general fiber is a disjoint union of P^s (cf. Corollary 4.2 (hi) below). 
Explicitly, say F4 = AxQi + A2Q2, deg Ai = deg Qi = 2, everything in sight as general as needed. 
Then 

F 4 = (A, - tQ 2 )Q 1 + {A 2 + iQi)Q 2 , 
so F4 contains the pencil of elliptic quartics {A\ — tQ 2 ,A 2 + tQi),t e P 1 ; setting t — 00, we find 
(Qi,Q2)- Similarly, get (Qi — tA 2 , Q 2 + tA\). This is one and the same pencil. But there is 
also (Ai — tA 2 ,Q 2 +tQi). In general, these 2 pencils are disjoint. Looking at them as curves in 
X = G(2,7 r 2 ), we actually get a Plucker-embedded conic, (Ai — tA 2 ) A (Q 2 + tQi) = A\ A Q 2 + 
t(Ai A Qi — A 2 A Q 2 ) — t 2 A 2 A Q\, disjoint from Y (see (3.1)). In particular, capping each conic 
against the Pliicker hyperplane class II = — c%A, we find 2. As before, we may write 

dcgNL(W,d)= J H 33 nNL{W 7 d). 

The cycle p 2 H 33 n NL(W, d) can be represented by a sum of deg NL(W, d) disjoint unions of pairs 
ofP x 's. Hence 

deg NL{W,d) = - [ U-H 33 nNL{W,d) = - [ U- ( Pl ) * H 33 DNL(W 7 d) 

= i/ s n.c I5(£l ). 

The latter integral can be computed via Bott's formula and we get 38475, cf. the script in [13]. 
This has been found independently in [3] with different techniques, using [10]. 

4. The fibers of p 2 
The main result needed to validate the above enumeration is the following. 

4.1. Proposition. Let C C P 3 be a smooth irreducible curve of degree e and genus g. Let (I > 
and let FcP 3 be a general surface of degree d containing C . Then C is the only effective divisor 
of degree e and arithmetic genus g on F . 

Proof. By [8, Cor. II. 3. 8] we have that Pic(F) is freely generated by its hyperplane section H and 
C. Let C be an effective divisor of degree e and arithmetic genus g on F. Then there are two 
integers a, b such that, on F, we have C ~ aH + bC. Now e = H ■ C = ad + be, so that 

(4.1) a= ^ 1 - & )" 
By adjunction formula we have 

(4.2) C 2 = 2 Pa (C) - 2 - K F ■ C = 2g - 2 - e(d - 4) = {C f < 
as d > 0. Now 

C 2 = (C) 2 = (aH + bCf = a 2 d + b 2 C 2 + 2abe 
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and using (4.1) we get 

(1 - b 2 )(e 2 - dC 2 ) = 0. 

Note that e 2 - dC 2 > by (4.2), whence b = ±1. If 6 = -1 we have from (4.1) a = ^ £ Z, as 
d > 0. Therefore we deduce that 6=1 and a = 0, that is C" ~ C. Since C 2 < we must have 

c = c. □ 

4.2. Corollary. Lei 6e an irreducible subvariety of a Hilbert scheme component of curves in P 3 
of degree e and arithmetic genus g with general member smooth. Let F N = P (H° (Op3 (d))) . Then 

(i) there is a do such that for all d > do the projection map p2 '■ NL(W, d) — > ¥ N is generically 
infective. 

(ii) // W is the family of elliptic quartics then we can take do = 5, i.e., p2 is generically one-to- 
one for d > 5 and 

(hi) for d = 4 ; the general fiber of NL(W,4) -— V NL(W, 4) is two disjoint P 1 's. 

Proof. We know from (2.4) that pi : NL(W,d) — > W is a projective bundle. Hence NL(W, d) = 
P2(NL(W,d)) is irreducible and a general element F G NL(W, d) can be identified with a general 
hypcrsurface of degree d containing a general C E W. Hence assertion (i) follows from the pro- 
postion. Assertion (ii) also follows, except for d = 8. In this case, with the notation as in (4.1), if 
b = — 1 we would get a = 1, whence C ~ H — C so that C would be contained in a plane, absurd. 
If b = 1, may proceed as at the end of the proof of the above proposition. 
If d = 4 then we get instead C ~ 2H — C. The exact sequence 

— >• O f — > C F (C) — )• 0c(C) = Oc — ► 

shows that \O f {C)\ = P 1 and similarly \O f (2H - C)\ ^ |Of(C")| = P 1 - Moreover there is no 
curve D on F such that D ~ C and Z? ~ 2iJ — C for then 2C ~ 2H giving the contradiction 
= C 2 = H 2 — 4. This proves that, in this case, the general fiber of P2 is two disjoint P^s, 

(4.3) P2\F) = \O f {C)\ U \O f (2H - C)\. 

□ 

4.3. The formula. In view of Prop. 2.1, it suffices to find the degrees of NL(W, d) for 3 • 16+1 
values of d > 5 and interpolate. This is done in [13]. We obtain 

( d 3 2 ) (l06984881d 29 - 34095 14775d 28 + 57226549167d 27 - 
643910429259d 26 + 5267988084411d 25 - 31628193518727d 24 + 126939490699539d 23 
-144650681793207d 22 - 2701978741671631d 21 + 28913126128882647d 20 - 
182919422241175163d 19 + 858473373993063183d 18 - 3061191057059772423d 17 
+7448109470245631187d 16 - 3841505361473930575d 15 - 80644842327962348733d 14 + 
568059231910087276234d 13 - 2560865812030993315212d 12 + 9159430737614259196104d n 
-27608527286339077691280d 10 + 71605637662357479581024d 9 
-160009170853633152594240d 8 + 303685692157317249665152d 7 
-473993548940769326728704d 6 + 571505502502703378479104d 5 
-459462480152611231457280d 4 + 111908571251948243582976d 3 
+251116612534424272896000d 2 - 328452832055501940326400d 

+136886449647246114816000^) j (2 27 ■ 3 9 ■ 5 2 ■ 7 2 • 11 • 13). 

4.4. Acknowledgment. Thanks are due to the referee for kindly pointing out to us a sizable 
portion of the argument that required clarification. We arc also grateful to Kristian Rancstad and 
Rahul Pandharipande for correcting an error in §3.2, and calling our attention to [3] and [10]. 
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APPENDIX 

Script for Singular. This builds on the original script for Maple found in P. Meurer, [11]. Copy 
and paste it into a Singular session. 

//<"frac4bott" ; 

option (noloadLib.noredef ine) ; 
LIB "ring. lib"; 

proc Proj (n){ 

executeO'ring P"+string(n)+"=0,x(0. .n) ,dp; 
def xx=ideal(x(0. .n)) ; export xx; 

keepring P"+string(n) ) ; nameof (basering) ; baser ing; }• 
proc omit(w,a) 

{list ll,lw,v;int i,j,s,t,ch; 
s, t=size( w),size( a); 

if(typeof( a)o"list") { 
list 11; 

if (typeof (a)=="poly") { 

t=l;ll[l]=a;> 

if (typeof (a)=="ideal") { 

t=size (matrix (a)) ;ch=l;} 

if (typeof (a)=="intvec" or typeof (a)=="int") 
{t=size ( (a) ) ; ch=l ; } 
if (ch==l){ 

for( i=l; i<= t; i++){ll[i]=a[i] ;»} 
else{ll= a;} 

if (typeof (w)<>" list" ){ 
list lw; 

if (typeof (w)=="poly") { 
s=l ; lw [1] =w; } 
if (typeof (w)=="ideal") { 
s=size (matrix (w) ) ;ch=l;} 

if (typeof (w)=="intvec" or typeof (w)=="int") 
{s=size ( (w) ) ; ch=l ; } 
if (ch==l){ 

for( i=l; i<= s; i++) {lw [i] =w [i] ; » 
def ty=typeof (w) ; 
}else{lw=w; } 

for(i=l;i<=s;i++) 

{ch=0; 
for(j=l; j<=t; j++){ 
if (lw [i] ==11 [j] ) {ch=l ; break ; » 

if (ch==0){ 
v=insert (v , lw [i] , size (v) ) ; } 

}if (def ined(ty)<>0){ 
if (size (v) <>0) {def st=string(v) ; } 
else{def st="0";} 
execute(ty+" vv="+st);} 



SURFACES CONTAINING AN ELLIPTIC QUARTIC CURVE 



9 



else{def vv=v;} 
return (vv) ; } 

proc origin(v,p) 

{def q=p ; if (typeof (q) <>"list " ) { 

for(int i=l ; i<=size (v) ;i++){q=subst(q,v[i] ,poly(0)) ;} 

return (q) ; } 

else{f or (int i=l ; i<=size (q) ; i++){q[i] =origin(v,q[i] ) ; } 
return (q) ; } > 

proc hilbp// (ideal I) 

{int nargs=size (#); ideal I=#[l]; 

//start comparing rings 

poly pvars=indets (I) ; 

if (pvars==0) {pvars=indets (vars (1 . . (nvars (basering) ) ) ) ; } 
if (pvarsOindets (vars (1 . . (nvars (basering) ) ) ) ) 
{if (nargs==3) 

{print ("too many vars in basering "+nameof (basering) + 
", will change to less");} 
if (nargs>=2){pvars=indets(#[2] ) ;} 
string st="ring r0=(0" ; if (npars (basering) <>0) {st=st+" , "+parstr (basering) ; } 
st=st+") , ("; 

for (int i=l ; i<=(size (pvars) -1) ; i++) {st=st+string(pvars [i] )+",";} 
st=st+string(pvars [i] )+" ) , (dp) " ; 
st=st+" ; ideal I="+string(I) ; 
execute (st) ; if (nargs==3) {where () ; }} 

intvec iv=hilbPoly (I) ; if (nargs==3) {print (iv) ;} 
int s=size (iv) ; string st.stl; 
for(int i=l ; i<=s ; i++){ 

if (iv [i] <>0){stl=string(iv [i] ) ; if (i==l){st=stl ; } 
else{if (i==2){if (iv [i] <>1) {stl=stl+"*t " ; } 

else{stl="t" ;}if (size (st) >0) {if (iv[i] <0) {st=st+stl ; }else{st=st+"+"+stl ; }} 
else{st=stl ; }} 

else{if (iv[i]ol){stl=stl+"*t~"+string(i-l) ;} 

else{stl="t~"+string(i-l) ;}if (size(st)>0){ 
if (iv [i] <0) {st=st+stl ; >else{st=st+"+"+st 1 ; » 

else{st=stl ;}»» 
if (s>=3){st=" (" + st + ")/" + string(f actorial(s-l) ) ;> 
return(st) ; } 

proc indets (p) //returns sum of variables 

{int i,j,s;poly v,w,vl;ideal pp; 

if (typeof (p)=="int" or typeof (p)==" intvec" 

or typeof (p)=="intmat" or typeof (p)==" number") 
{v=0 ; return ( v) ; } 
else{ 

if (typeof (p)=="poly") { 

for (i=l ; i<=nvars (basering) ; i++) 
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{if (reduce (p , std(var (i) ) )<>p) 
{v=v+var (i) ; } 

} 

return (v) ; } 

else{if (typeof (p)=="ideal") { 

pp=simplif y (p , 2) ;s=size(pp) ;w=0; 
if (s==0){v=0;return(v) ;} 
else{ 

for(i=l;i<=s;i++) 

{v=indets (pp [i] ) ; 
if (v<>0) {vl=leadmonom(v) ; 
w=w+v*vl~i ; }} ; 
return(indets (w) ) ;}}- 
else 

{if (typeof (p)=="vector" or typeof (p)=="matrix") 
{return (indets (ideal (p) ) ) ; 
} 

else 

{if (typeof (p)=="list") 

{def s=size(p);if (s==0){return(poly(0)) ;} 
else{ideal pp; 

f or (i=l ; i<=s ; i++) {pp [i] =indets(p [i] ) ; 

}return(indets (pp) ) ;}} 

}}}» 

proc vars(v) 

{ideal z;def vv=v;if (typeof (vv)=="intvec") { 
for (int i=l ; i<=size (vv) ; i++){z=z+var (vv [i] ) ; } 
return(z) ; } 

else{if (typeof (vv)=="ring") { 

if (string(vv)Ostring(basering) ) { 

def goback=nameof (basering) ; 

setring vv;def p=vars (intvec (1 . . nvars(vv) ) ) ; 

execute ("setring "+goback) ;def p=imap(vv,p) ; 
return (p) ; }else{return(vars (intvec (1 . .nvars(vv) ) ) ) ; } 
>» 

proc mdc 

{int s , i ; s=size (#); string ty=typeof (# [1] ) ; 
if (s==l){ 

if (ty=="poly" or ty=="int" ) {return(# [1] ) ; } 
else{ 

if (ty=="list" or ty=="vector" or 
ty=="ideal" or ty=="intvec" 
or ty=="matrix" ) 

{if (ty=="vector" or ty=="matrix") 
{ideal ii=ideal (# [1] ) ; 
} 
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else{ 

execute(ty+" ii=#[l];"); 

} 

if (ty<>"intvec"){ 
ii=simplify(ii,2) ;} 
i=size (ii) ; 

if (i==0){return(poly(0)) ;} 

else{return(mdc(ii [1 . . i] ) ) ; 

} 

> 

} 

> 

if (s==2){ 

return (gcd(#[l] ,#[2])) ; 

} 

if (s>2){ 

execute (ty+" p=#[l];"); 
p=gcd(#[l] ,#[2]); 
for(i=3;i<=s;i++){p=gcd(p,#[i] ) ;} 
return (p) ;}} 

proc dotprod( wl , w2) { 
if (typeof (wl)=="list" 
& typeof (wl)=="list") 
{int s=min(size (wl) ,size(w2) ) ;poly p; 
for (int j=l; j<=s; j++){p=p+wl[j]*w2[j] ;} 
return (p) ; } 
else{list vl,v2; 

if (typeof (wl)=="vector") {//"//wl vector"; 
for (int i=l ; i<=nrows(wl) ; i++){vl [i] =wl [i] ;}} 

else{ if (typeof (wl)=="ideal") {//"//wl ideal"; 
for(int i=l;i<=ncols(wl) ; i++) {vl [i] =wl [i] ;» 

else{if (typeof (wl)=="intvec"){//"//wl intvec" ; 
for(int i=l;i<=size(wl) ; i++) {vl [i] =wl [i] ; }} 

else{ERROR(typeof (wl)+ 
"?wl?not (int)vector nor ideal");}}} 
if (typeof (w2)=="vector") {//"//w2 vector"; 
for (int i=l;i<=nrows(w2) ; i++) {v2 [i] =w2 [i] ;}} 

else{ if (typeof (w2)==" ideal") {//"//w2 ideal"; 
for (int i=l;i<=ncols(w2) ; i++) {v2 [i] =w2 [i] ;}} 

else{if (typeof (w2) ==" intvec " ) {// " //w2 intvec " ; 
for (int i=l;i<=size(w2) ; i++) {v2 [i] =w2 [i] ;}} 

else{ERROR (typeof (w2)+ 
"?w2?not (int)vector nor ideal");}}} 
return (dotprod (vl ,v2) ) ;}} 

proc min//intvec or list or whatever? 
{def s=size (#) ; 
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def m=#[l] ;for(int i=2 ; i<=s ; i++) {if (# [i] <m) {m=# [i] ; » 
return (m) ; } 

proc prod// (ideal y) 

{ def s ,ty=size (#) ,typeof (# [1] ) ; 

poly p=l;int i; 

if (ty==" ideal") {def y=#[l]; 

matrix m=matrix(y) ; //bypass size problem 

f or (i=l ; i<=ncols (m) ; i++) {p=p*m[l , i] ; } 

} 

else{ 

for(i=l;i<=s;i++){p=p*#[i] ; }} 
return (p) ; } 

proc coefmon(poly p, poly m) 
//coeff of monom m in poly p 
{ if(m==l or p==0){return(jet(p,0)) ;} 

else{ideal ml=std(m) ; def v=std(pol2id(indets (m) ) ) ;poly qO; 

for(int i=l ; i<=size (p) ; i++){def q=reduce (p [i] ,ml) ; 

if (q==0){q=p[i]/m; 

if (reduce (q,v)==q) {qO=qO+q; }}} 
return (qO) ; } } 

proc pol2id(poly p){ 
ideal I;for(int i=l ; i<=size (p) ; i++) 
{I[i]=p[i] ;} return(I) ;} 

proc sumf rac// (list ab,list cd,...) 
{def nargs=size (#) ; //"//"+string(nargs) ; 
if (nargs==2) {list ab=#[l];list cd=# [2] ; 
def 1= (ideal(ab[l]*cd[2]+ab[2]*cd[l] , 
ab[2]*cd[2]));def p=gcd(l[l] ,1[2]); 
returndist (1 [1] /p, 1 [2] /p) ) ; } 
else{def ll=sumf rac (# [1] ,#[2] ) ; 
for(int i=3;i<=nargs;i++){ 
ll=sumfrac(ll,#[i] ) ;> return(ll) ; » 

proc subf rac (list ab.list cd) 
{def 1= (ideal(ab[l]*cd[2]-ab[2]*cd[l] , 
ab[2]*cd[2]));def p=gcd(l[l] ,1[2]); 
returndist (1 [1] /p, 1 [2] /p) ) ; } 

proc mulf rac (list ab,list cd) 

{def 1= ideal (ab[l]*cd[l] ,ab [2] *cd[2] ) ; 

def p=gcd(l[l] ,1[2]); 

returndist (1 [1] /p, 1 [2] /p) ) ; } 



proc invf rac (list ab) 
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{return(list(ab[2] ,ab[l] )) ;} 

proc divf rac (list ab.list cd) 
{return (mulfrac (ab , invf rac (cd) ) ) ; } 

proc mon(vs,d) 

{if (typeof (vs)=="poly"){def vvs=pol2id(vs) ;> 
else{def vvs=vs;} 
return(vvs'd) ;} 

proc mylead(poly h){ 
def oldring=nameof (basering) ; 
def goback="setring "+oldring; 
def hl=h; 

if (f ind(string(hl) , "x")<>0){setring rx;> 
else{setring rz;def xx=zs;} 
executeC'def f =imap("+oldring+" ,hl) ") ; 
def iv,a=leadexp(f ) ,leadcoef (f ) ; 
a=int (a) ; 

f =(dotprod(iv,xx) ) ; 

def al ,hl=string(a) ,string(f ) ; 

if (f ind(hl , "x")==0) { setring rxz; 

executeC'def f="+hl); 

f=z2x(f ) ;hl=string(f ) ; setring rx; 

executeC'def f="+hl); 

} 

f =x2p(f ) ;list mm; 

executeC'int j = "+al); 

for(int i=l;i<=(j) ;i++){mm[i]=f ;} 

if (oldring<>"rx"){execute(goback) ; 

def mm=imap(rx,mm) ;} 

return (mm) ; } 

proc pesos (H){ 

//H= poly or list(num,den) of fraction 

list rres; 

def ty=typeof (H) ; 

if (ty=="poly" or ty=="number" or ty=="int"){ 
if (indets (H) ==0) {return (list (0) ) ; }else{ 
rres=mylead(H [1] ) ; 

for(int i=2; i<=size(H) ; i++){ 
rres=rres+mylead(H [i] ) ;} 
return(rres) ;}} 

else{if (ty=="list" and size(H)==2) { 

def num,den=H[l. .2] ; 

//den assumed monomial 

rres=pesos (H [1] ) ; def denn=pesos (den) ; 

for(int i=l ; i<=size (rres) ; i++) { 
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rres [i] =rres [i] -denn[l] ;} 
return(rres) ; } 
else{if (ty==" ideal") 
{rres=pesos (H [1] ) ; 

for(int i=2; i<=size (matrix(H) ) ; i++) { 
rres=rres+pesos (H [i] ) ; }} return(rres) ; }}■} 

proc cherns(H){ 

if (def ined(t)o-l){def t=var(l) ;} 
poly sigma= 1; 

for(int j=l ; j<=size(H) ; j++){ 

sigma=reduce ( (1+H [j] *t) *sigma, std(t~ (DIM+1) ) ) ; } 
list 1; 

for (int j=l ; j<=DIM; j++){ 
1 [j] =coefmon( sigma, t~j ) ; } 
return(l) ; } 

proc topchern(H) {return (prod (H) ); } 

proc TGrass(idd_, E_){ 

if (typeof (idd_)o"list") { 

def En_,iddn_=list(E_,l) , list (idd_ , 1) ;} 

else{def iddn_,En_= idd_,E_;} 

def p=subf rac (En_ , iddn_) ; 

def L=divf rac (p , list (iddn_ [1] [1] , iddn_ [2] ) ) ; 
for(int j=2; j<=size(iddn_[l] ) ; j++){ 

L=sumf rac (L.divf rac (p, list (iddn_[l] [j] ,iddn_[2] ))) ;} 
return (L) ; } 

proc x2p(f){ 

if (def ined(n)==0) {return ("//need int n");} 

if (def ined(ivP)==0){ 

return ("//need intvec ivP=0 ,2 ,7 , 10 . . . ") ; } 

if (def ined (x (n) ) ==0) {return ( 
"//need Proj (n) ; rx=Pn; ideal xx=x(0 . .n) ") ; } 
def g=subst(f ,x(0) ,poly(ivP [1] ) ) ; 
for(int i=l ; i<=n; i++){ 
g=subst(g,x(i) ,poly(ivP[i+l] )) ;> 
return (g); } 
//fixpts start here 
//<"4ellsing" ; <"f rac4bottarXiv" ; 

int n=3; Proj (n); def rx=P3;ring r=0,t,dp;r=r+rx;setring r; 
imapall (P3) ; 

def xx2=xx~2;def xx2s=sum(xx2) ; 

def xx4=std(xx2~2) ; def xx4s=sum(xx4) ; 

def xs=sum(xx) ; 
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intvec ivP=0, 1, 5, 18;int DIM=16; 
list G; 

for(int i=l ; i<=size (xx2) ; i++) { 
for(int j=i+l ; j<=size(xx2) ; j++){ 
G[l+size(G)]=list(list(xx2[i] ,xx2 [j] , 
"G:ql,q2,tg"), 
TGrass(xx2[i]+xx2 [j] ,xx2s) 

);}> 

list Z,G2; //Z=P3* x G(2,4), plane, line 
for(int i=l;i<= size (G) ; i++) { 
def l=G[i]; def ql ,q2=l [1] [1] , 1 [1] [2] ; 
def p=gcd(ql ,q2) ; 

if( pOl ) {def ll,12=ql/p,q2/p; def 
tg=sumf rac (TGrass (11+12, xs) ,TGrass(p,xs)) ; 
def nml=subf rac (1 [2] ,tg) ; 
Z[l+size(Z)] 

=list (list (p, 11, 12, i, "Z:p=mdc, 11,12, pos.G.tg.nml" ) , 
tg.nml) ;} 

else{def q4=ideal(ql,q2)*xx2; 
G2[l+size(G2)] 

=list(list(ql,q2) ,G[i] [2] ,q4) ; 

if( size(q4)<>19 ){ string(i) ;break; >}} 

list E1,G2E1,W; //exc over Z 

//should get welldef system of 8 cubics everywhere 
//W=P2-bdle over clsd orbit of Z 
for(int i=l;i<= size (Z) ; i++) { 

def l=Z[i]; def L1=G [1 [1] [4] ] [1] ;L1=L1 [1 . . 2] ; 

nml=l [3] ; 

def num,den=nml [1 . . 2] ; 

for(int il=l;il<= size (num) ; il++) { 

def exc=sumf rac (list (0 , 1) , list (num [il] ,den)) ; 

int jl=0; 

for(int j=l ; j<=size(Ll) ; j++){ def ll=Ll[j]; 

def j0=ll/exc[2] ; 

if(j0<>0 Hjl=l; 

//check div instead 

break; }} 

if (jl==l){ 

def L2=omit(Ll ,11) ; 

L2=L2[1] ,ll+t*exc[l]*jO; 

def qs=sat (ideal (L2[l. . 2])*xx,t) [1] ; 

qs=reduce (qs , xx4) ; 

qs=std(origin(t,qs)) ;def j0=size(qs) ,hilbp(qs) ; 

tg=divfrac(subfrac(nml,exc) ,exc) ;//(nml-exc)/exc 
tg=sumf rac (tg, 1 [2] ) ; //tg f iber+base 
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tg=sumfrac(tg,exc) ;def 11=1 [1]; 

El [l+size(ED] =list(list(ll[l. .4] ,i, 
"El: mdc,ll,12,pos.G,pos.Z,2:tg,3:exc,4:8 cubics"), 
tg, exc.qs); 

if (j0[2]o"4*t") { def L=mdc(qs); //L=l[l][l] 
def qsO=ideal(qs/L) ; def 11=1 [1] [2] +1 [1] [3] -L; 
def tgO=sumf rac (TGrass (L ,xs) //tg P3dual 
,TGrass(ll,xs-L)) ; //tg fiber WO— >P3dual 
def qO=origin(ideal(L,ll) , list (sum(qsO) ,xx2s)) ; 
tgO=sumf rac (tgO, TGrass (qO[l] ,q0[2] )) ; 

W[l+size(W)]=list (list (L,ll ,i, size (El) , 
"W:mdc,ll,pos.Z,pos.El,2:tgW,3:nml,4:exc,5:8th q,6:8 qs"), 
tgO, subfrac(tg,tgO) , exc , qO [1] , qsO) ; } 
else {def qs4=std(qs*xx) ; 
G2El[l+size(G2El)] 
=list(list(size(El) ,tg,qs4)) ; 
if( size(qs4)<>19 ){ "//?"+string(i) ;} 
» 

else//( jl==0 ) 

{ "//??"+string(i , " , " , il) ; }//never happens, no need 
>} 

def h = "4*t"; 

list E2; //exc over W; //P8-bdle over W; dimW=3+2+2 

for(int i=l;i<= size(W);i++) { 

def l=W[i] ; 1=1 [1] ; 

def L,ll,q0=l[l. .2] ,W[i] [5]; 

def nml=W[i] [3] ; 

def num , den=nml [1 . . 2] ; 
def qs=El[l[4]] [4] ; 
qs=std(qs*xx) ; 

for (int il=l;il<= size (num) ; il++) { 

exc=sumf rac (list (0, 1) , list (num [il] ,den) ) ; 

tg=divfrac(subfrac(nml,exc) ,exc) ;//(nml-exc)/exc 

tg=sumfrac(tg,W[i] [2]);//tg fibra+base 

tg=sumf rac (tg , exc) ; 

def qsl= qs + num[il]; 
E2 [l+size(E2)]=list(list(L,ll,qO,i, 
"E2: L,ll,qO,pos.W,2:tg,3:exc,4:19 qtics"), 
tg, exc, qsl); 
» 

def l=list(G,Z,El,E2,G2,G2El) ; 

setring rx; 

def l=imap(r,l); 

def G,Z,E1,E2,G2,G2E1=1[1. .size(l)] ; 



//actual numerical computations at fixpts start here 
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//case d=4 

poly f ;def X=G2; 

for (int i=l;i<=size(X) ;i++) { 

def H=-sum(pesos (sum(X [i] [1] ) ) ) ; //1st chern class 
//=Plucker hypln class 

def qs , c=std(X [i] [3] ) ,topchern(pesos (X [i] [2] ) ) ; 
def d=cherns (pesos (kbase (qs , 4) ) ) [DIM-1] *H ; 
f=f+d/c; 

if(i mod 10 ==1) {"//in G2, remain "+string(size(X)-i) ; 

system("sh","date"); 

» 

"//done G2"; 

X=G2E1;//X[1] [1] . . . 

for(int i=l;i<=size(X) ;i++){ 

def l=El[X[i] [1] [1]] [1] ; 

def ql2=l[l]*l[2]+l[l]*l[3] ; 

def H=-sum(pesos (ql2) ) ;//lst chern class 

def qs,c=std(X[i] [1] [3] ) .topchern (pesos (X[i] [1] [2] ) ) ; 

def d=cherns (pesos (kbase (qs , 4) ) ) [DIM-1] *H ; 

f=f+d/c; 

if(i mod 10 ==1) {"//in G2E1 , remain "+string(size(X)-i) ; 
system ("sh" , "date") ; } } 

X=E2; 

for(int i=l;i<=size(X) ;i++){ 
def l=X[i] ; 

def ql2=l[l] [1]~2+1[1] [1]*1[1] [2] ; 

def H=-sum(pesos (ql2) ) ;//lst chern class 

def qs , c=std(l [4] ) , topchern (pesos (1 [2] ) ) ; 

def d=cherns (pesos (kbase (qs , 4) ) ) [DIM-1] *H ; 

f=f+d/c; 

if(i mod 10 ==1) {"//in E2, remain "+string(size(X)-i) ; 

systemC'sh" , "date") ; } } f; 

f=f/4; 

"//case d=4: "+string(f ) ; 

//choose maxdeg>=5 ; need 48... 
int minideg,maxideg=54,54; 

def X=G2;poly f (minideg. .maxideg) ; 
for(int i=l;i<=size(X) ;i++){ 

def qs , c=std(X [i] [3] ) ,topchern(pesos (X [i] [2] ) ) ; 
for (int il=minideg; il<=maxideg; il++){ 
def d=cherns (pesos (kbase (qs,il))) [DIM] ; 
f (il)=f (il)+d/c;> 
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if(i mod 10 ==1) {"//in G2, remain "+string(size(X)-i) ; 
system("sh","date"); } }"//done G2"; 

X=G2E1;//X[1] [1] . . . 

for(int i=l;i<=size(X) ;i++){ 

def qs,c=std(X[i] [1] [3] ) .topchern (pesos (X[i] [1] [2] )) ; 
for (int il=minideg; il<=maxideg; il++){ 
def d=cherns (pesos (kbase (qs , il) ) ) [DIM] ; 
f (il)=f (il)+d/c;> 

if(i mod 10 ==1) {"//in G2E1 , remain "+string(size (X) -i) ; 
system("sh","date"); } }"//done G2E1"; 

X=E2; 

for (int i=l;i<=size(X) ;i++){ 

def qs,c=std(X[i] [4]) ,topchern(pesos (X [i] [2] ) ) ; 
for (int il=minideg; il<=maxideg; il++){ 
def d=cherns (pesos (kbase (qs , il) ) ) [DIM] ; 
f (il)=f (il)+d/c;> 

if(i mod 10 ==1) {"//in E2, remain "+string(size(X)-i) ; 
systemO'sh" , "date") ; } }"//done E2"; 

for (int il=minideg; il<=maxideg; il++){ 
string (il,",",! (il));} 

if (l>2){//check answer 

def t=var (1) ; 

def 

rr= (t-2) * (t-3) * (t-4) * (106984881*t-29-3409514775*t ~28+ 

57226549167*t-27-643910429259*t-26+5267988084411*t-25 

-31628193518727*t-24+126939490699539*t-23-144650681793207*t~22 

-2701978741671631*t-21+28913126128882647*t-20 

-182919422241175163*t-19+858473373993063183*t-18-\ 

3061191057059772423*t~17+7448109470245631187*t~16 

-3841505361473930575*t-15-80644842327962348733*t'14+ 

568059231910087276234*t-13-2560865812030993315212*t-12+ 

9159430737614259196104*t-ll-27608527286339077691280*t-10+ 

71605637662357479581024*t-9-160009170853633152594240*t~8+ 

303685692157317249665152*t~7-473993548940769326728704*t~6+ 

571505502502703378479104*t~5-459462480152611231457280*t~4+ 

111908571251948243582976*t-3+251116612534424272896000*t-2-\ 

328452832055501940326400*t+136886449647246114816000) 

/(number (2) "28* (3) "10* (5) ~2* (7) "2*11*13) ; 

//subs(t=54,rr) ; 

}1>2 
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